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ABSTRACT 

We analyze the spatial distributions of young stars in Taurus- Auriga and Upper Sco as determined 
from the two-point correlation function (i.e. the mean surface density of neighbors). The corre- 
sponding power-law fits allow us to determine the fractal dimensions of each association's spatial 
distribution, measure the stellar velocity dispersions, and distinguish between the bound binary pop- 
ulation and chance alignments of members. We find that the fractal dimension of Taurus is D ~1.05, 
consistent with its filamentary structure. The fractal dimension of Upper Sco may be even shallower 
[D ~0.7), but this fit is uncertain due to the limited area and possible spatially- variable incomplete- 
ness. We also find that random stellar motions have erased all primordial structure on scales of <0.07° 
in Taurus and ^1.7° in Upper Sco; given ages of ^1 Myr and ^^5 Myr, the corresponding internal 
velocity dispersions are ~0.2 km s~^ and ~1.0 km s~^, respectively. Finally, we find that binaries can 
be distinguished from chance alignments at separations of <120" (17000 AU) in Taurus and <75" 
(11000 AU) in Upper Sco. The binary populations in these associations that we previously studied, 
spanning separations of 3-30", is dominated by binary systems. However, the few lowest-mass pairs 
(Mprim^O.3 Mq) might be chance alignments. 

Subject headings: stars:formation, stars:kinematics, stars :pre~main sequence, methods:statistical 



1. INTRODUCTION 

The spatial distribution of young stars is a powerful 
diagnostic of their formation and early evolution. Young 
stars trace the gas distribution from which they formed, 
so the large-scale structure of a young association retains 
these primordial features after the gas has been accreted 
or dispersed. On intermediate scales, the absence of 
structure indicates the typical distance over which stars 
have randomly dispersed since their birth, and therefore 
the velocity dispersion for the association. Finally, the 
enhanced stellar density on small scales outlines the bi- 
nary population, distinguishing bound binary systems 
from chance alignments between young stars. Some of 
these topics have been addressed in previous work on 
young star distributions (Gomez et al. 1993; Larson 
1995; Simon 1997; Bate et al. 1998; Hartmann 2002), 
but the modern census of several key star-forming re- 
gions is more complete and extends to lower masses than 
a decade ago, so the analysis is worth revisiting. 

The traditional tool for studying spatial distributions 
is the two-point correlation function (hereafter TPCF). 
The TPCF, 'w{9), is defined as the number of excess 
pairs of objects with a given separation over the ex- 
pected number for a random distribution (Peebles 1980). 
The TPCF is proportional to the mean surface density 
of neighbors, so it is often recast in terms of this more 
intuitive quantity: S(6') = {N^/A)[l + w{9)], where A is 
the survey area and A^* is the total number of stars. 

In this letter, we describe an updated relation for Yj{6) 
in Taurus and present the first such analysis for Up- 
per Sco, then we fit power laws for the different angu- 
lar regimes. Finally, we interpret our results to address 
three questions: What is the primordial fractal dimen- 
sion of star-forming regions, and how does it relate to 
their observed geometry? What is the primordial veloc- 
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ity dispersion suggested by each association's randomiza- 
tion? And what is a wide binary companion, and can it 
be distinguished from an unbound association member? 

2. THE CORRELATION FUNCTIONS OF TAURUS AND 
UPPER SCO 

We compiled our Taurus sample from the member list 
in Kraus & Hillenbrand (2007a, 2008; hereafter KH07a 
and KH08), plus the Class O/I sources that were com- 
piled by Kenyon & Hartmann (1995). We omitted the 
latter sources from our multiplicity surveys because their 
stellar properties are uncertain, but we include them here 
because that information is not necessary for clustering 
analysis. We have also included the partial list of new 
sources identified in data from the Taurus Spitzer Legacy 
Project (Padgett et al. 2006) as described by Luhman et 
al. (2006). For separations of <30", we have calculated 
the surface density of neighbors only among those sources 
included in our initial wide binary survey. We have ne- 
glected the Class O/I and heavily embedded sources be- 
cause only some have been surveyed for multiplicity in 
the mid-infrared wavelengths (e.g. Duchene et al. 2004), 
and not with uniform sensitivity. Our full sample consists 
of 272 members, while the binary-regime sample consists 
of 226. The Taurus sample is almost certainly incom- 
plete, as a number of additional candidates have been 
identified in the Taurus Legacy Project (Padgett et al. 
2006) and the XEST survey (Scelsi et al. 2007). How- 
ever, preliminary reports suggest an increment of ^20% 
in the total sample. Even if these new members do not 
trace the known distribution, their influence should be 
modest. 

The census of Upper Sco across the full association 
is very incomplete, so we implemented our analysis for 
intermediate and large separations {0 > 30") using 
only members in two heavily-studied fields observed by 
Preibisch et al. (2002), the 2dfE and 2dfW fields. The 
census of members in these fields is not complete, but we 
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TABLE 1 
Power Law Fits 



Regime Sep Range'' a C (log deg ^) 







Upper Sco 






Binary 


3"-30" 




-1.44±0.41 


2.98±0.12 (at 9.5") 


Intermediate 


2.8'-1.5° 




-0.12±0.02 


1.537±0.010 (at 16') 


Association 


1.5°-4.7° 




-1.31±0.09 


1.174±0.011 (at 2.6°) 






Taurus 




Binary 


3"-30" 




-1.53±0.32 


3.28±0.10 (at 9.5") 
1.62±0.05 (at 2.8') 


Intermediate 


1.6'-5.0' 




0.12±0.39 


Association 


5.0'-4.7° 




-0.951±0.007 


0.650±0.005 (at 1.2°) 



There is a small range of separations between the binary and intermediate regimes where the data are consistent with our power law fits, 
but the uncertainties are too large for those data to contribute meaningfully to the fits. 
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Fig. 1. — Locations of stars in Taurus and Upper Sco, superimposed on 60/im IRAS images. Members are denoted by green crosses, 
while the sample fields in Upper Sco are denoted by blue circles. The field of view is 17° in Taurus and 5° in Upper Sco. Known members 
in Upper Sco outline the dusty clouds in the northern field, suggesting systematic incompleteness for extincted members. 



expect that it is the least incomplete. As for Taurus, we 
calculated the surface density of neighbors at <30" us- 
ing the full sample of our wide binary surveys; this choice 
maximizes our sample size for small separations (where 
the statistics are weakest). The 2dfE/2dfW and binary 
samples consist of 162 and 352 members. 

In Figure 1, we plot the locations of our sample mem- 
bers superimposed on archival 60//m IRAS images. In 
Upper Sco, we see evidence of incompleteness for the 
northern field. Most of the known members outline the 
dusty regions, suggesting that any members in these re- 
gions were too extincted to have been identified. As we 
discuss later, this could affect the TPCF on scales of 
In Taurus, the distribution traces the filamentary dust, 
though there are also many filaments that do not include 
any known members. 

We directly measured S(6') for Taurus because our 
sample spans the entire area of the association. However, 
for bounded subsets (as in Upper Sco) , it is often easier to 
evaluate the TPCF via a Monte Carlo-based definition, 
w(Q) = Np(e)/Nr{e) - 1, where Np{e) is the number of 
pairs with separations in a bin centered on 9 and Nr{0) 



is the expected number of pairs for a random distribu- 
tion of objects over the bounded area (Hewett 1982). 
The advantage is that this method does not require edge 
corrections, unlike direct measurement of S(0). In both 
cases, we report our results as S(^) since it is a more visu- 
ally motivated quantity than wlO). In Figure 2, we plot 
I](0) for Upper Sco (top) and Taurus (bottom) spanning 
a separation range of 3" to 10°. 

Based on the predicted time-evolution of young associ- 
ations (Bate et al. 1998), we expect that S(0) can be fit 
with a twice-broken power law, corresponding to struc- 
ture on three scales. At small scales, bound binary sys- 
tems yield a relatively steep power law. At large scales 
(and for young ages, <1 crossing time), intra-association 
clustering yields a shallower (but nonzero) power law 
that corresponds to the primordial structure of the asso- 
ciation. Finally, at intermediate separations, the random 
motion of association members acts to smooth out the 
primordial structure and yield a constant surface density 
(and thus a slope near zero, according to the simulations 
of Bate et al. 1998). The first knee (transition between 
gravitationally bound multiplicity and a smooth random- 
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Fig. 2. — Two-point correlation functions for members of Upper Sco and Taurus. These plots show the surface density of neighbors 
as a function of separation, with B in degrees (bottom axis) or in parsecs (top axis). The observations are from our recent wide 

binary survey (KH08; filled circles) or membership surveys in the literature (open circles). For each association, we have fit power laws 
to the small-scale regime (red; binary systems), the large-scale regime (blue; association members distributed according to the primordial 
structure), and the intermediate regime (green; association members with a randomized spatial distribution). 



ized distribution) corresponds to the maximum angular 
scale for distinguishing binary systems, while the second 
knee (transition between a random distribution and pri- 
mordial structure) corresponds to an angular scale that 
depends on the age since members were released from 
their natal gas clouds, r, and the internal velocity disper- 
sion, Uini, where 9(XTVint- Hartmann (2002) suggested 
that this break also could indicate the mean spacing of 
cores along filaments (the Jeans length), which assumes 
that stars have randomized by a smaller angular scale 
and that the inferred value characteristic angular scale, 
the inferred value of Vint is an upper limit. 

In Table 1, we summarize our weighted least-squares 
fits for the power law slope a and zero-point C in each 
regime. The binary regime was fit in the range probed 
in our survey of wide multiplicity (3-30"), while the in- 
termediate and association regimes were fit in the ranges 



where the error bars were ^3%. We established the zero 
point of each fit at the logarithmic center of the angular 
range in order to minimize correlation between CTq and 
<Tc- In Upper Sco, both the inner and middle power laws 
are clearly defined, but the fit for the outer regime is un- 
certain because the angular scale is similar to the size of 
the survey area (~2-4°). In Taurus, the inner and outer 
power laws are clearly defined, but the fit for the inter- 
mediate regime is uncertain. The TPCF at separations 
of 2-4' is flat and diverges from the fit for larger and 
smaller separations by 3-5tT, so we provisionally assume 
that this separation range represents the intermediate 
regime. The points at smaller separations also fall be- 
low the projection of the association-regime power law, 
while the points at larger separations agree well with the 
overall fit, suggesting that our inferred value of Vmt is 
at most an upper limit. The locations of the first knee. 
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where the two power hiws axe equal, are (?i.[75co^75"and 
di,Tau^^'20"; the respective locations of the second knee 
are ^2,!7Sco~l-7° and 02,Tau~O.O7°. The formal uncer- 
tainties in these measurements are only ~2-3%, but the 
errors are dominated by systematic uncertainties in the 
membership census and in the angular range over which 
to fit each regime. 

3. ASSOCIATION REGIME: THE FRACTAL DIMENSION 
OF TAURUS 

The primordial spatial distribution of young stars 
shoiild trace the overdensities in the original gas distri- 
bution from which those stars formed (e.g. Hartmann 
2002; Bate et al. 2003). Even if these gas distributions 
have dispersed, the remnants of primordial structure in 
the stellar distribution can still provide a key constraint 
to the distribution of overdensities during star formation. 
Early studies of TPCFs have suggested that current (and 
presumably primordial) stellar distributions are fractal in 
nature (e.g. Larson 1995: Simon 1997), with self-similar 
structure on a range of angular scales. Similar TPCFs 
can be reproduced (at least over a decade of separation) 
with simpler distributions like a finite number of non- 
fractal subclusters following a simple profile (Bate 
et al. 1998). However, our TPCF for Taurus follows 
a single power law across >2 decades of separation, so 
it appears to be genuinely self-similar. The dimension- 
ality D oi a fractal distribution indicates the extent to 
which it fills space, such that the number of neighbors A'^ 
within a distance 9 goes as N{6)(x9^ . This parameter is 
related to the surface density of neighbors; if E(0)oc0", 
then D = a+2 (Larson 1995). 

The fractal dimension is a result of the turbulent 
fragmentation that leads to star formation, and most 
models yield filamentuary structure (i.e. a dimension 
near unity). As we showed in Section 2, the observed 
power-law slope for Taurus in the large-scale regime is 
a= — 0.951±0.007, indicating that the fractal dimension 
{D = 1.049±0.007) is indeed close to unity. This result 
is consistent with visual inspection of the stellar distri- 
bution, as well as with CO maps of the remaining gas 
distribution (e.g. Goldsmith et al. 2008). Our value is 
significantly lower than the fractal dimension suggested 
by Larson (1995) and Simon (1997), D = 1.4, but close 
to the more recent value suggested by Hartmann et al. 
(2002). Our sample is significantly more complete than 
the older samples; based on our reconstruction of those 
samples, most of the new (typically low-mass) members 
are located near the major concentrations (e.g. Strom 
& Strom 1994; Briceno et al. 2002) rather than in the 
more distributed population (e.g. Slesnick et al. 2006). 
These members increase the surface density of neighbors 
at small separations, yielding a steeper slope for T,{9). 
However, the census is still incomplete (Section 2) and if 
the incompleteness is spatially variable, such as for heav- 
ily embedded brown dwarfs, then our updated power-law 
slope could be incorrect. 

We are hesitant to estimate the fractal dimension in 
Upper Seo. The appropriate regime in the TPCF in- 
cludes only two separation bins, so the choice of bin lo- 
cations could significantly affect the slope. Incomplete- 
ness in the dusty northern region could also influence the 
inferred large-scale structure. However, if we adopt our 
power-law fit from Section 2 (a= — 1.31±0.09), we find 



that D = 0.69±0.09 on scales of -2°. 

4. INTERMEDIATE REGIME: THE PRIMORDIAL 
VELOCITY DISPERSION 

The angular scales over which structure has been ran- 
domized, as indicated by the location of the second knee 
in T,{0), directly constrains the primordial velocity dis- 
persion for each association (e.g. Bate et al. 1998). This 
constraint is particularly important for low-density asso- 
ciations like Taurus and Upper Sco because the expected 
velocity dispersion (<l-2 km s^^; Frink et al. 1997) may 
be too low to be measured easily via a direct method (like 
high-resolution spectroscopy to determine radial veloci- 
ties). OB and T associations are not bound once their 
unaccreted gas is expelled (e.g. Lada et al. 1984), so the 
internal velocity dispersion is critical for determining how 
long they can persist as recognizable moving groups (like 
the (3 Pic, TW Hya, or /x Oph associations; Webb et al. 
1999; Zuckerman et al. 2004; Mamajek 2006) and how 
long substructure can remain in these moving groups. 

Allowing for projection effects, the angular scales of 
each TPCF's outer knee correspond to physical disper- 
sion scales of ~0.23 pc in Taurus and ~6 pc in Upper Sco. 
Given the characteristic ages of each association (^1 Myr 
and ~5 Myr, respectively), the corresponding character- 
istic velocity dispersions are ~0.2 km s~^ and ~1.0 km 
s""'^. As we previously discussed, there is uncertainty in 
the fits, so this values should be taken with caution. We 
also note that these values represent the velocity disper- 
sion with respect to other stars only within an angular 
distance of ~^fe„ee- We can't rule out the possibility that 
larger substructures are moving coherently with a higher 
velocity dispersion, only that any substructure with an- 
gular size 9 is not moving with sufficient speed {9^9 /t) 
that its angular displacement from birth is of order 0. 
This limit also suggests an explanation for the larger ve- 
locity dispersion in Upper Sco; even if the velocity dis- 
persion within ~0.1-0.2 pc substructures is the same as 
in Taurus, the observed TPCF eoidd be reproduced if the 
velocity dispersion between those substructures is ~1 km 
s~^. A scale dependence in the velocity dispersion could 
also explain previous proper motion studies in Taurus, 
which found velocity dispersions within the major sub- 
clumps (on scales of ~l-3 pc; e.g. Jones & Herbig 1979) 
that were ~1 km s~^. 

A similar effect has been noted in locations like the 
ONC, where radial velocities show an overall north-south 
gradient of ~5 km s^^ in addition to the local velocity 
dispersion of 2-3 km s^^ (Furesz et al. 2008). However, 
there is also observational evidence that small-scale ve- 
locity dispersions are higher in denser clusters; submil- 
limeter observations of IRSl in NGC 2264 (Wilhams & 
Garland 2002) found that six protostellar cores (span- 
ning 0.44 pc) had a velocity dispersion of 0.9 km s~^, 
which much higher than the velocity dispersion that we 
find in Taurus, though also closer to the value for scales 
of 1-3 pc suggested by Jones & Herbig (1979). 

Our results suggest that regions like Taurus and Up- 
per Seo are even less dynamically active, relative to the 
ONC, than their lower densities might imply. The veloc- 
ity dispersions also provide a direct estimate of the virial 
velocity in the natal environment (before the removal 
of gas) and therefore jointly constrain the typical mass 
and size of a star-forming clump: Mu^~(3GM^)/(5ii) 



5 



or M/i?-(5'i;2)/(3G), yielding M/R -15 in Taurus and 
M/R ~550 in Upper Sco, where the mass is in solar 
masses and the radius is in parsecs. Thus, the primor- 
dial star-forming structures that are now dispersing with 
these characteristic velocities were smaller and/or more 
massive in Upper Sco than in Taurus. 

5. BINARY REGIME: WHAT IS A BINARY SYSTEM? 

The existence and properties of wide binary systems 
are critical for constraining multiple star formation in the 
limiting case of large separations and early times. If wide 
binaries form out of a single protostellar clump, then the 
maximum separation also constrains the maximum size 
of clumps that can collapse to become bound systems. As 
previous authors have suggested (e.g. Larson 1995), the 
outer edge of the young binary separation distribution 
is similar to the mean Jeans length for nearby molecu- 
lar clouds. This limit is also similar to the maximum 
separation seen in the field (e.g. Duquennoy & Mayor 
1991), suggesting that some wide binaries join the field 
without being subjected to significant dynamical inter- 
actions. However, study of young binaries is complicated 
by the difficulty of distinguishing gravitationally bound 
binary pairs from coeval, comoving association members 
that are aligned in projection. We addressed this issue 
for a single system in Upper Sco (UScoJ1606-1935; Kraus 
& Hillenbrand 2007b) by calculating the association's 
TPCF to determine the probability that it is a bound 
system; we now extend our analysis to the full known 
populations of Taurus and Upper Sco. 

We find that the transition between the binary and in- 
termediate regimes occurs at -^11000 AU in Upper Sco 
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and ~17000 AU in Taurus. The binary population there- 
fore extends at least to these angular scales, but we can 
not distinguish binary companions from chance align- 
ments outside this limit. The difference between these 
regimes is a result of the higher total wide binary fre- 
quency in Taurus (KHOTa, KH08), as the overall sur- 
face density of "contaminant" co-association members is 
similar in both associations. The number statistics do 
not support any assertions regarding the outer maximum 
limit of binary formation, but this angular scale matches 
both the maximum binary separation seen in the field 
and the typical Jeans length, so we do not expect to find 
many binary systems with wider separations. 

Candidate companions inside this limit could also be 
chance alignments, but the probability drops for progres- 
sively smaller separations. In Upper Sco, we expect —3.3 
chance alignments with separations of 15-30" from an 
intermediate- or high-mass member {Mprim >0.4 Mq), 
plus another —2.4 chance alignments of two low-mass 
members. The number of high-mass chance alignments 
is far lower than the total number of pairs, which sug- 
gests that the vast majority are bound binaries. How- 
ever, our wide binary survey found only four wide pairs 
of low-mass companions, so it is unclear whether any 
are genuine binary systems. The contamination rate is 
moderately lower in Taurus, yielding <1 contaminant in 
either mass range, but the results are similar. Most of 
the high-mass pairs are binary systems, but the two low- 
mass pairs may or may not be bound binary systems. 
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